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00 ! Abstract 

G^ ■ 

— , In this paper a new formalism based on exterior differential systems 

^ ' is derived for perfect-fluid spacetimes endowed with an abelian orthogo- 

I I nally transitive G2 group of motions acting on spacelike surfaces. This 

5-H ■ formulation allows simplifications of Einstein equations and it can be ap- 

plied for different purposes. As an example a singularity-free metric is 

^ , rederived in this framework. A sufficient condition for a diagonal metric 

•'^ ■ to be geodesically complete is also provided. 

b : 

- - ■ 1 Introduction 

Perfect-fluid spacetimes endowed with an abelian G2 group of isometrics have 
been used for describing many different physical situations. When the group 
is acting on timelike surfaces, they have been extensively used for describing 
axisymmetric compact objects in stationary rotation (cfr. for a review). 

On the other hand, if the group acts on spacelike surfaces, the applications 
are different ||]. A classification of these spacetimes is given in Q. They can 
model spacetimes where two plane gravitational waves are colliding (cfr. |^ 
for a review) , but they are also useful for describing inhomogeneous cosmologies 
(cfr. 1^,1^ for a review) in an attempt of coping with the inhomogeneity present 
in our Universe. An interesting feature of G2 cosmologies is the possibility of 
avoiding initial and final singularities (cfr. , Q , for a review) and therefore 
physics will be valid in the whole spacetime. These models satisfy the causality 
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and energy conditions and just fail to contain trapped sets, according to the 
well-known singularity theorems ||^ , . 

In this paper we shall try to cope with two features concerning nonsingular 
perfect fluid orthogonally transitive G2 cosmological models: First of all one 
has to devise a method for obtanining exact models and then one should check 
whether the solution is singular or not. Both of these will be the aim of this 
paper. 

A new tetrad formalism based on differential forms will be introduced for 
deriving results about spatially inhomogeneous spacetimes. It will be shown that 
the methods initially devised for stationary axially symmetric spacetimes 
are also useful when the group of motions acts on spacelike surfaces. The 1-forms 
that are used in this formalism will be shown to have kinematical meaning when 
considering spatial congruences. If the congruence corresponds to an invariantly 
defined quantity (for instance, in the cylindrical case, the axial Killing), the 
kinematical properties can be used to classify the solutions. Also, using the 
remaining gauge freedom, the tetrad can be adapted to the congruence so that 
the exterior system can be written in terms of 1-forms with an invariant and 
physical interpretation. Moreover, when one tries to obtain an exact model, it is 
generally useful to impose certain assumptions on these kinematical quantities 
(now with an invariant interpretation). As the approach is grounded on an 
exterior differential system, the coordinates can be chosen according to the 
Ansatz that is performed, instead of fixing them from the beginning. 

Assuming one has obtained a new cosmological model, it usually takes 
lengthy calculations to determine if there are singularitities in it. The ques- 
tion is settled if the curvature invariants are already singular, but if they are 
regular in the whole spacetime, there is a priori no reason to assume that every 
geodesic is complete. In fact, there are cosmologies with regular curvature in- 
variants that are incomplete and therefore singular, despite no quantity becomes 
unboundedly large Therefore we deem convenient to have a general result 
that may simplify the task of analyzing the issue of geodesic completeness. This 
matter will be addressed for the diagonal case in this paper. 

Let us describe in more detail the contents of the paper: 

In the second section spacelike congruences in a general spacetime will be 
studied in order to achieve an interpretation for the mathematical quantities 
that will appear in the formalism. In the stationary axisymmetric case timelike 
congruences were considered and their tangent field could be considered as a 
velocity and therefore the interpretation was straightforward. Another way of 
interpreting spacelike congruences is due to Greenberg JT^ , but in this paper a 
different approach will be followed. 

In the third section the formalism is written in terms of an exterior system 
of equations that include Cartan and Einstein field equations as well as their 
integrability conditions. The set of equations will be simplified taking advantage 
of the remaining gauge freedom. 

As an example of how the exterior system can be used for obtaining exact 
solutions the singularity-free model in will be obtained within the formalism 
in the fourth section. 
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The question of geodesic completeness of diagonal inhomogeneous cosmolog- 
ical models will be addressed in the fifth section and a theorem will be derived 
as a sufHcient condition for a model to be nonsingular. 

This condition will be shown to be weak enough to comprise all known di- 
agonal singularity-free inhomogeneous cosniological models in the sixth section. 



2 Spacelike congruences 

The kinematical properties of a timelike congruence can be defined by decompos- 
ing the covariant derivative of the timelike vector field defined by the congruence 
[ l3 , ||l5| . An analysis of the spacelike congruences has been done by Greenberg 
[12|. In his analysis Greenberg introduces an observer moving with a four ve- 
locity w'^. Projecting the covariant derivative of the vector field defined by the 
congruence orthogonal to w'^ he obtains "spacelike" quantities characterizing 
the congruence. In this section we will follow a different approach to study the 
kinematical properties of a spacelike congruence. We will use a straightforward 
translation of the timelike congruence analysis. In this way we obtain kinemat- 
ical properties not necessarily with a spacelike character, but we will find that 
these kinematical quantities can be used in a natural way to formulate a sim- 
plified differential form approach for spacetimes with two commuting spacelike 
Killing vectors. 

Let us consider a congruence of spacelike curves, 

a;"=x"(y^T) , a =1,2, 3, 

where the three constants = specify a particular curve and r is the arc- 
length. We can define a unit tangent vector, 

n'^ = , rian = 1, 

dr 

We will assume in the following that the congruence defines, at least locally, 
a vector field n^{x°'). 

Let us take a particular curve (C) in the congruence, specified by three 
constants (y"), and a point p in the curve (C) characterized by an arc-length r. 
Now consider another curve in the congruence (C*) near (C) and specified by 
constants y"* — y°' + 5y°' and a point p* in (C* ) with the same arc- length (r) 
as p. Then, up to first order, we have, 

which, in general, is not orthogonal to n". In order to get a vector orthogonal 
to we introduce the projector tensor, 
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such that Pl^nf^ = 0. Then we can define, 

6^x" = F/^fe^ {6±x°'nc, = 0). 

Note that S±x°' can be spacehke, timeUke, or null. 

The rate of change of the connecting vector of two spacelike curves of the 
congruences allows us to characterize locally the congruence. First, let us con- 
sider the change in the 'modulus'. Assume S^x" is timelike or spacelike, then, 

6(61)^ ^ g^pS^x'^S^x^ , e = ±l. 

The rate of change of 61 along the congruence is, 

{6iy ^ S^x'^S^x^ 1 

+ (1) 



where we define. 



SI "' SI (5/ 3 ' 



$ ^ (2) 



n„ = rin-fin 



c^iin" (3) 

Sap = "-(q;/3) - "(q^I/?) - g^^'a/S, (4) 



where S* is a trace-free symmetric tensor. 
If we define the unitary connecting vector, 



SI 

the equation (nl) can be written as follows. 



= eSape'^e^ + i$. (5) 



SI 3 
On the other hand for the rate of change of e" we get, 

P'3je")- = {W^^ -f ~ e^^.e^e'-Oe" (6) 

where, 

Wap = n^a-p] - n[anpY 

If e" is an eigenvector of S^^^ (S'^'^e" = Ae^) then the previous equations 
reduce to, 

P^Je")-= 14^/^6". (7) 

It is interesting to note that if we know tIq, Sap, Wap, and $ in a given 
point we can reconstruct the congruence locally using that, 

na-p = riaUp + Wap + SaP + ^^Pa/S- 
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Some important properties are, 

P^pn'^ = , n„n" = , Wo^pn^ = = So^pn^ . 

In order to have a better characterization of the spacelike congruence we 
can study the eigenvalue problem for the trace- free symmetric three-dimensional 
tensor SajS that can be formulated as follows, 

{S^p - XPo,p)v^ - 0. 

As the quadratic forms defined by Sap and Pap are not definite forms, a 
standard analysis gives us the following different situations: 

• Ai, A2, and A3 = — Ai — A2 are three different real numbers. In this case 
we have one timelike eigenvector and two spacelike eigenvectors. They are 
mutually orthogonal. 

• Ai, Ai are complex conjugate and A2 = — 2Rc(Ai). In this case we have 
two complex conjugate eigenvectors, m and m corresponding to the com- 
plex eigenvalues and a real spacelike eigenvector orthogonal to m and to. 
The complex eigenvector can be written as rn = a + ib, where a is timelike 
and b spacelike and they are mutually orthogonal and normalized to —1/2 
and 1/2 respectively. 

• Ai = A2 ^ 0, A3 = — 2Ai are real numbers. In this case there is a 
null eigenvector corresponding to the double eigenvalue and one spacelike 
eigenvector orthogonal to the plane containing the null eigenvector. 

• Ai = A2 = A3 = 0. In this case there is a null eigenvector corresponding 
to the triple eigenvalue. 

The interpretation of the kinematical properties of a spacelike congruence 
depends on the physical interpretation of the vector field n and the character 
of the eigenvalues an eigenvectors of Sap. Assume for instance that the three 
eingenvalues are real and different. As a consequence the three eigenvectors are 
mutually orthogonal and there is one timelike and two spacelike eigenvectors. 
In order to give an interpretation in this case, let us consider a curve (C) in the 
congruence and a point (P) in it corresponding to r = tq. Assume that in the 
spacelike subspace in the orthogonal hyperspace to (C) in (P) we have a disk. 
Orthogonal to it there is a unit timelike vector. If we imagine different curves 
of the congruences crossing the border of the disk and the final point of the 
timelike vector, going from tq to tq + St, the disk changes to an ellipse and the 
direction and modulus of the timelike vector also changes. These changes are 
determined by the kinematical variables of the congruences <&, Sap, and Wap 
and the final effect can be decomposed in the following two steps: 
Step 1: ($ and Sap effects) 

The influence of (f> and Sap on the disk and timelike vector can be understood 
using equation If e" is an eigenvector of Sap and verifies Sape^ — Ae" the 
equation reduces to, 

{SIX . 1 . 



Hence, the two directions indicated by the two spacehke eigenvectors will trans- 
form to the principal directions of the ellipse and the length of the principal axis 
will be obtained from the previous equation. The new direction of the timelike 
vector will be given by the direction of the timelike eigenvector and the modulus 
will be obtained from the previous equation. In this way we can think of the 
spacelike congruence as generating a "virtual tube flux" where the transverse 
section is given by the ellipse, the velocity of the fluid by the timelike unit eigen- 
vector of Sai3 and its modulus can be interpreted as the density, whose changes 
are given by the previous equation (note that, in principle, the virtual fluid has 
nothing to do with any real perfect fluid. The physical interpretation of such 
virtual fluid depends of the physical interpretation of the congruence) . 
Step 2: {Wa(3 effects) 

Finally the effect of Wap on the "tube flux" is to "bend and twist the tube" . 
This result follows from the equation (j^). 

The only remaining kinematical variable is fia but this represents the cur- 
vature of the selected curve in the congruence. 

In the case when Saj3 has degenerate or complex eigenvalues the interpre- 
tation is not so straightforward. We will study those cases in the particular 
situation when the metric admits two Killing vectors. 

2.1 Kinematical Variables in a tetrad formalism 

Using a tetrad adapted to n such that 6^ = n, 

ds^ = -6*" + + 9"^ ^e"^ + ® 0^, 

the covariant derivative of n takes the following form, 

na-b = -7206, 

where 7" be = ~Oi-j^b^c (where {ca] is the orthonormal frame dual to {9°"}) are 
the Ricci rotation coefficients. The kinematical properties of the congruence of 
n read as follows. 



n 



= -7020 + 7121 - 7233 (9) 

Wt = (7021 - 7120)^" A + (7023 + 723o)f° A + (7123 + 723i)f ^ A 9^ (10) 

2 11 12 1 

St = (g7020 + g7i2i - 27233)6*° 6*° + (g7020 + g7i2i + 27233)6*^ fXi 6*^ + 

112 

(37020 - 37121 - 37233)6!^ ®9^+ (11) 

(7021 + 7120)6*° «)s 9^ -f (7023 - 723o)6l° ®s 0^ + (7123 - 7231))6l^ ®s 6'^ 

where ®s is the symmetric part of the tensor product, that is, v — {v ® 

w + w ® v)/2. 
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3 Tetrad formalism for a spacetime endowed with 
two spacelike orthogonally transitive commut- 
ing Killing vectors 

If a spacetime has two spacelike Killing vectors 77} defining an abelian or- 
thogonally transitive G2 group of motions we can choose the tetrad for the 
spacetime {0° ,9^ ,9^ ,9^}, 

ds^ ^ -~9" ®9^ + 9^®9^+9^®9^ + 9^® 9^, 

such that 6*2 and 9^ are in lm{^,r]} and 0. 

L^6^ = = Lf^O^ ; -^c^^ = = LijO^ 1 

L^9^ = = L,,6|2 ; L^9^ = = £^6*^ 

It is important to notice that there is a residual 50(1, 1) gauge in the {0" — 9^} 
subspace and also a 5*0(2) gauge in the {9^ — 9^} subspace. 

For the tetrad presented above a family of independent non-vanishing Ricci 
rotation coefficients are the following, 



7010, 7011: 7022, 7023 — 7032,7033,7122,7123 — 7l32 , 7l33, 7230, 7231 ■ 



Then, the kinematical variables reduce to, 

fl = 7022^° + 71220' (12) 

$ = (13) 

Wt = (7o23 +723o)f" A 6l3-H (7123+7231)6'' A 6*3 (14) 

St = (7023 - 723o)f°'^s 6*^ + (7123- 7231) )f'®sf^ (15) 
It is interesting to define 1-forms a, lo, and cr, 

a = -70226*" - 71226'' (16) 

w = -(7023 + 7230)6*" - (7123 + 7231)6'' (17) 

= (7230 - 7023)6*" + (7231 - 7123)6'' (18) 

such that the kinematical variables can be written as. 



h = —a 
Wt = -LoA 9^ 
St = ~<J (»s 9^ 



so that a, w, and a completely parametrize the kinematical variables. 

The eigenvalue problem for St simplifies and there are three different cases. 
Defining A = 7023 - 7230 and B = 7123 - 7231 we get: 
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1. B = eA where = 1 (cr null) 

In this case ct is a null form a ~ A{d'^ + e6^). We have Ai = A2 = A3 = 
whose eigenvector is proportional to a. 

2. B > A {a spacelike) 

In this case we have Ai=0, Vi (x *a (timelike). 
A2 = y/W^ , V2(xa + - ^26*3 (spacelike). 
A3 = -\/S2 - , ^3 cx CT - V-B2 - ^26*3 (spacelike). 

3. B < A {a timelike) 

In this case we have Ai = , Vi oc *cr (spacelike). 

A2 = WA^ ~ ^2 ^ V2 0ca + WA^~B^e^. 

A3 = ~WA^ - 52 , 1/3 CX (T - iV^2 _ 52g)3 

In order to complete the family of 1-forms that will be used to write an 
exterior system equivalent to the Einstein equations, we introduce two new 1- 
forms (3 and v, 

P = -(7022 + 7033)^*° - (7122 + 7133)^*^ (19) 
= 7oioe° + 7oiie'. (20) 

The interpretation of these 1-forms will be found immediately below. 
3.1 Vanishing torsion equations 

With our choice of tetrad and variables, the vanishing torsion equations can be 
written as. 



de° = (21) 

dd^ = vhe° (22) 

de^ = ahO'^+uhO^ (23) 

dO^ = (/3 - a) A 61^ + (T A (24) 



The meaning of v and (3 is now clear, since u is the connection form in the 
e° - subspace, while d(6'2 A 9^) ^ (3 A 9"^ A 9^ shows that /3 describes the 
expansion of the volume in the 9^ — 9'^ subspace []. 

^Note that /3 is proportional to the differential of the transitivity surface area element, as 
^fined in and also is proportional to the differential of the function W in eq. (15.3) in 
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3.2 First Bianchi identities 



The first Bianchi identities are the integrability conditions for the equations 
( plfp^ ). By exterior differentiation of those equations we find, 

dp = 0, (25) 

dn + KAS = 0, (26) 

d(5 - K A r2 0, (27) 

dK. + nAS = 0, (28) 

where fl = P — 2a, 6 = uj + a and k = uj — a. 
3.3 Einstein field equations 

We consider the energy-momentum tensor of a perfect fluid, T^^ = (fi+p)uaU^ + 
pgau (where Ua is the velocity of the fluid, /i is the energy density and p the 
pressure of the fluid), as the source of the gravitational field. In appropriate 
units, the Einstein equations read as, 

or, equivalently, 

Rau = {n+p)UaUy + - {pl.-p)ga„. 

After writing the Ricci tensor in the orthonormal coframe in terms of the 
Ricci rotation coefficients and identifying the kinematical 1-forms and their 
derivatives, the Einstein field equations can be combined to produce the follow- 
ing equivalent exterior system,^ 

d*n + PA*n + KA*6 = 0, (29) 
d*5 + P A*5- kA*VL = 0, (30) 
d*P + PA*P = {^i-p)uA*u, (31) 

di/ - A */3 + A *0 + i(5 A *(5 = -^(^ + p)u A (32) 
d*P+]^PA*P+'UlA*h + ]^6A*5 + 2vAP = -{ji+p)u A^u, (33) 
dp +]^P AP +^^1 AVL + ]^6 A~5 + 2v A^P = -{ji+p)uAu, (34) 

where the tilde operation on 1-forms in the 6^ — 9^ subspace is a reflection 
with respect to the direction determined by f?^. If A = XqO'^ + Ai^^, then 
A = \q9^ — \i6^ . The * operator is the Hodge dual in the 2-subspace 9^ — 9^ . 

similar calculation can be found in for the stationary axisymmetric case. 
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3.4 Integrability conditions 



As u and the other physically relevant 1-forms have vanishing Lie derivative with 
respect to the two Killing fields, the coefficients of the forms we have considered 
are independent of the variables in the 6^ — 9^ subspace. The integrability 
conditions for the system of equations that we have considered above are either 
trivial or already incorporated in the system, except for the contracted Bianchi 
identity, 

Tf"" = 0, 

which in our case reduces to, 

du-\ —dpAu = 0, (35) 

d*u+(l3-\ —diJ,)A*u = 0, (36) 

and therefore the fluid is irrotational. 

We also have to take into account the following constraint indicating that 
the velocity of the fluid is a unit timelike vector: 

uA*u = -e° A0\ 



3.5 Simplification of the equations 

Recall that there is a remaining freedom of choice consisting in a rotation of 6^ 
and 9^ in their own plane, as well as a rotation of 0'^ and 9'^ in their own plane. 
This freedom can be used to simplify the differential system introduced above. 

3.5.1 Gauge freedom in the 9^ — 9^ subspace 

Under a rotation through an angle if, 




a2 



■sm{ip) cos((/?) ) \ 9^ 
the kinematical 1-forms transform as. 



n \ _ / cos(2(^) -sin(2(^) \ f ^ 
6 ) ^ [ sin(2(^) cos(2</?) ) [ S 



(37) 



P = f3, (38) 
k = K + 2dip, (39) 



(40) 



Without loss of generality, we can use this freedom in order to impose, for 
instance, 

k = eS ; (e^ = 1), 
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(note that k = S -i^ a = and k = —5 cj = 0) or, 

fi = en ; (e^ = 1) 

as the integrabihty condition for those equations derived from eq. (|3^) are iden- 
tically satisfied using the first Bianchi identities. A more general gauge fixing 
condition is, 

k — cos(2iy9o)^^ — sm(2ipo)5 ; ipQ ~ const., 

and in the same way as in the previous cases the integrabihty conditions are 
satisfied using eqs. ( p5[]2^ ). 

It is important, in order to characterize the solutions of the equations, to 
know if a given solution is "diagonal" or not. In an adapted gauge this is 
equivalent to S = = k. In a non-adapted gauge we get dn — and A (5 = 0. 
In other words if we have a solution satisfying the previous equations, with an 
adequate transformation, the metric can be written in a diagonal form. 

3.5.2 Gauge freedom in the 6'^ — 9^ subspace 

The freedom in the 9^ — 9^ subspace can be used to align 6*" with a timelike 
kinematical form or 9^ with a spacelike kinematical form; a good candidate 
seems to be u but a careful analysis shows that /3 is more useful. In this way 
we have to separate three different cases: 
1) (3 null 

If /3 is a null form then we have, 

*/3 = e/3 ; = 1. 

Using this relation in eq. ( |33| ) and adding eq. ( |33| ) with eq. (|4|) we get 

nA{n + e*tt)+SA{S + e*S)+ Apu A (?i + e * it) = 0, 
and introducing the relations, 

n = na9° + ni9\ 

S = So9" + di9\ 

in the previous equation we get, 

(f7o + e^if + {5o + e5if + 4p(wo -1- tuif = 0. 

As a consequence if p > we get that f2, (5, and u are null forms. Hence, 
if (3 null then either p < or u is null. 

The cases with p < or u null are not very physical for a perfect fluid and 
we will not consider them here. The only remaining case is p = but from 
( ^l| ) and the fact that /3 is null we get ^ — and then the solution represents a 
vacuum spacetime. 
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2) (3 timelike 

If f3 is timelike we can align 6*° with it. With this choice, the equations in 
[ pl}|36| to be solved reduce to the system composed by: 

u/\^u = e^^(3 A */3, 

plus the last two vanishing torsion equations, (|23| ) an d (p^ , the first Bianchi 
identities, [2^^ , the first three Einstein equations, |29| - t3l| , the intcgrability 



conditions, (|35| ) and (36), and, finally 



3 1 

dQ = -I3~ -e^^{{ji-p)l3+{ji + p)[< (3,u> u+<(3,*u> *m] (41) 

+^[<(3,n> n+<f3,*n> *n] + ^[<(3,s> s+<f3,*s> *s]}, 

where <, > stands for the scalar product and < 6'", 6'" —1, < 9^, 9^ >— 1 , 
and < 9^, 6^ >= =< 9^, 0^ >. It is important to notice that the connection 
can be solved algebraically from (^) and (^J) and the i' thus obtained satisfies 
( ^2[ ) identically, using the remaining equations of the system. Thus, we can 
forget about v, which may be obtained trivially after we have solved the system 
presented above. A useful expression for v is the following, 

iy = *dQ-[l- {n-p)e^'^] * p. 

Also, it has to be noticed that for a generic 1-form a in the 9° — 6'^-subspace 
we can write, 

a = — < a,u > u+ < a, *u > *u, 

and as a consequence, 

o^Bu =< a,u > u+ < a,*u > *u, 

can be interpreted as a reflection of a with respect to u. An analogous inter- 
pretation is possible for the similar expression changing u by SI (or S ) . Finally 
let us note that the intcgrability condition for the equation (|4^) is identically 
satisfied using the rest of the system. 
3) (3 spacelike 

If (3 is spacelike we can align 9^ with it. In a similar way as in the timelike 
case, the system can be reduced to [E3[«^, pBj36| and, 



(3 = e-«0i 

3 1 

dQ = -/?+ -e2'3{(^-p)/3+ (^ + p)[< /3,u > u+</3,*M> *m] (42) 

+^[<(3,n> n+<i3,*n> *n] + ^[< i3,s > s+<f3,*s> *s]}. 
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The reasoning for timclike /? is valid for the spacehke case. The 1-form v 
can be written as, 

v^^dQ-\\ + (/i - p)e2<3] * (5. 

Again the integrabihty condition for the equation (^) is identically satisfied. 
Components of the differential of Q in the light-cone and positivity 

The spacelike or timelike character of /3 is essential in many senses |2], pTf . 
Here we present another one related with the gradient of the modulus of /3, 
parametrized by e^'^ . 

Hence, let us study the components of the differential of Q in the light-cone, 
that is, along the null geodesic directions on the 0" — 0^-subspace: 
/3 timelike 

We have, 

/3 = e-^e" 

S = 6a9° + Si9^ 
u = UQe^ + uie"" 

dQ = doQe" + diQe^ 



and then we get, 

P spacelike 

The only new equation is, 

and then the new components are, 

i(f}o±17i)2 + l(5o±^i)2 

We can see that the right-hand side of the equation for {di±do)Q is positive 
in the spacelike case (when the energy condition /i > \p\ is satisfied) and has not 
definite sign in the timelike case. Therefore in the spacelike case the transitivity 
cylinders are not trapped surfaces. 
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4 An example 



In order to show how the formahsm works let us derive the nonsingular solution 
in [^Sj. This nonsingular cosmology is diagonal, therefore we can choose a 
gauge such that uj = a — 0. The 1-form /? is spacelike as in every other regular 
inhomogeneous spacetinie and 9^ is chosen parallel to it as indicated in section 
3.5.2. 

Concerning the fluid, two Ansatze are made. The velocity of the fluid is taken 
to be —6*" and the fluid is stiff, that is, p — ^. Then we get that u = (3. 
Introducing this expression in (p5|) and (Bfl) the resulting equations imply, 



dg = -i(iln^ + /3. (43) 

The equations are written using a 1-form basis formed by the closed forms, 
a and (3. Locally this means that we can write, 

a = du, P = dv, (44) 

and therefore ( ^3| ) can be integrated to yield, 

/,e2« = aV^ (45) 

where a is a constant. 



Now we have to solve ( 29 ) and ( pOD (note that the only remaining equation 
( p2[ ) can be integrated by line integrals, when we have solved the rest of the 
equations, as its integrability condition is identically satisfied using the rest of 
the equations of the system) . Then, we introduce two functions /, g of the 
chosen variables it, v to express the Hodge dual of f3 as, 

du + f dv , , 

= , (46) 

9 

that allows us to complete the dual of our 1-form local basis, 

^^^^-fd^^Hl^f^^ (47) 
9 

taking into account the properties of the Hodge dual. 

Now (|2^) and ( ^ ) are equivalent to the following two partial differential 
equations, 

df fdg dg 
au ou ov 

_ ^ + _ ^ 0, (49) 
dv du du ' 

A simple solution of the system of equations can be obtained by requiring 

that / be a function of v. Then g can be integrated from (|49|), 



k-2ul (50) 
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where ft, is a function of v. Introducing this expression for g in (Q) we get, 

which can be easily integrated to get, 

f=P + qe'\ (52) 

h = -2q^e^'' -Ae^^'pqv + wqe'"', (53) 
where p, q, and w are constants. This allows us to write down the function (7, 



9 



= \/-2(j2e4'' -4:e^'"pqv + wqe^'' -Aque'^''. (54) 



The last equation to be integrated is (^). Using the expression for Q that 
has been obtained previously, it can be written as follows, 

^ = (2/2 - 2g^ - 2a^e'''')dv + (2 + Af)du, (55) 

which allows us to obtain the following expression for the fluid energy density, 

fi = k exp (^2p'^v + ^ q'^e'^" ^Ae^'"pqv - wqe^"" - a^e^" + 2u + Apu + Aque'^'" 

(56) 

where the constant k will be taken to be a? without loss of generality. 

The metric can be written in a simple form if we perform a change of variables 
and the coordinate v is taken to be Inr. If t is chosen accordingly so that the 
coordinates are isotropic, 

= Inr, u = — — (/r^ — plnr — (7^2, (57) 

the fluid density takes the form. 



H = a^ ^-2p(p+i) ^„ 1 ^2^4 _ (^2 _^ q + 2 pq)r^ -2qf - Apqf - 4 qh'^r'^^ , 

(58) 

after eliminating unnecessary parameters. 

The metric functions or the fluid density are singular unless p = 0,-1. Both 
cases are indeed the same solution just changing q for —q. In p3| , q is named 
13 and a? is a. 

Note that, in principle, this metric could also have been obtained by the 
generation algorithm due to Wainwright, Ince and Marsh [^ . 
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5 Geodesic Completeness 



In this section we shall attempt to give a sufhcient condition for a diagonal 
cylindrical G2-cosmology to be non-singular. Once one encounters a cosmologi- 
cal model with regular curvature invariants it usually takes lengthy and tedious 
calculations to prove that the metric is causally geodesically complete. The con- 
dition wc state here shall be large enough to comprise most known cylindrical 
non-singular cosmological models. 

For diagonal cylindrical G2-cosmology we can choose a gauge where lo = 
= (T (note that there is a residual gauge, a rotation with Lp — const.). The 
equations (|2^) and (|6|) can be integrated obtaining, 

a = -df , ^- (59) 

Using these expressions for a and /3, as well as the fact that w = = cr, we 
can integrate equations (|2^) and ( p4[ ) to get, 

0^ = e-^'dz, 
= pef d(j). 

All known non-singular diagonal cylindrical G2-cosmologies have a spacelike 
p. Hence, following the method described in section 3.5.2, we align 9^ with /3, 

(?i-e«^. (60) 
P 

As the cosmological fluid is irrotational, the velocity of the fluid can be 
written as, 

u — e^dt, 

using isotropic coordinates t and r such that, 

* dt = —dr and * dr = —dt, (61) 
that provides the following expression for 9''^ ^ — * 9^, 

= e'^{^dt+^dr). (62) 
P P 

The constraint u /\ *u = —e^^ [3 A imposes a relation between g and Q, 

e^Q = (63) 

pI - pi 

Using the previous results the metric takes the following form, 

ds^ = e^sit.r) |_^^2 + rf^2| ^ p2(^^^)g2/(t,r)^^2 ^ g-2 /(t,r)^^2^ (g^) 

Hence, the calculations will be performed in a chart where t and r are 
isotropic coordinates for the subspace orthogonal to the Killing orbits and </) 
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and z are coordinates adapted to the Killing fields. The usual ranges for these 
coordinates are chosen, 



-oo<i, z<oo, 0<r<oo, O<0<2 7r. (65) 

We shall assume from the beginning that the metric functions f,g,p are 
and that p is positive. Certainly a requirement is needed in order to avoid 
a singular Riemann tensor, but for geodesic completeness would be enough 
since just the affine connection is involved in the equations. 

Following 1^ , if we are to have a regular symmetry axis on the locus where 
^ — = 0; then we have to impose that, 

hm g(g''^dA, gradA) ^ 

P^o 4A ^ ' 

on approaching the axis. In our case this means, if the metric functions /, g are 
at the axis. 



p— ►O 

From now on wc shall denote the derivatives with respect to r and t by 
subscripts. 

There is no loss of generality in imposing that the axis is located on r = 0, 
since we are always free to choose a different set of isotropic coordinates T, R 
by performing a transformation, 

Tt = Tr = Ru (68) 

which amounts to take a solution of the one-dimensional wave equation, 

Rtt - Rrr = 0. (69) 

with a boundary condition i? = on the axis. This problem is underdetermined, 
since no initial condition has been imposed. For instance, one could take R ~ r 
and Rt ^ at t = 0. 

In order to determine whether the metric is geodesically complete, we shall 
write the expressions for the geodesic equations, 

F + T)kX^x'' = 0, (70) 

where the dot denotes derivative with respect to the affine parameter of the 
geodesic. 

In principle we would have to write a second order system of four differential 
equations, but two of them become first order due to the existence of isometrics. 
If u = {t, r, (j), z) is the 4- velocity of the geodesic and ^ is a Killing field then 
p = ^ • u is a constant of geodesic motion. The following quantities are then 
conserved along geodesies, 
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(71) 



p = ^-2f(t,r)^^ (72) 

respectively the angular momentum and the z component of the linear momen- 
tum of a unit mass test particle in free fall. 

There is also a conserved quantity, S, which takes the value zero for lightlike, 
one for timelikc and minus one for spacelike geodesies. This quantity arises from 
the fact that the 4-velocity u is normalized when an affine parametrization is 
used. We shall consider just future causal geodesies, 

The remaining second order equations read, 



i + gt{t, r)P + 2 grit, r)ir + gt{t, r)f^ - P'e^ {/(*.-)-s(*.'-)}/,(t, r) + 

-2{/(t,r)+s(t,r)} 



r + ff,(i,r)i2 +25t(i,r)tr + g,(<,r)f2 +PV^/(*''-)-s(*''-)>/,(t,r) - 

-2{f{t.r)+g{t,r)} 

-L^ {Pr{t, r) + Pit, r)frit, r)} , (75) 

pit, ry 

after substituting the derivatives of the cyclic coordinates 2, ^ for their expres- 
sions in terms of the constants of geodesic motion. 

These equations can be written in a more compact form making use of the 
equation for S, 



„-2g{t.r) 
|g2s(t,r)^|. _ ^ J p2g(t,r) 



5 + p2g2/(t,r)^^2, 



j-2/(t,r) 

7(vy 



=0, 



(76) 



e-2s(t,r) 



3-2/(t,r) 



= 0, (77) 



An important family of geodesies are lightlike radial geodesies (L = P 
i5 = 0), for which there is a constant k such that. 



i=\r\=ke-^<'^*'^\ 



(78) 



The derivative of the time coordinate must not grow too fast if these geodesies 
are complete. A sufficient condition is achieved by imposing that t does not grow 
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faster than exponentially for large values of t. This amounts to the following 
condition on g, 



9{t,r) > --\n\t + a\+b, 



(79) 



where a and b are constants. If this condition is fulfilled, t is defined for arbi- 
trarily large values of the affine parameter. 

The system of second order equations (fz^), { fn\ ) is shown to be equivalent 
to a system of first order equations by introducing a new function ^. Since, 



r „-2/(t,r) 
^2 _ .2 ^ e-2g(t.r) g ^ L^i— — + p2g2/(t.r) 

it is tempting to parametrize t, r by means of hyperbolic functions of f , 



(80) 



/ p-2/(t,r) 

t = e-a(t-'-\ 6 + + P^e^f(t,r) cosh^(i,r), 

V P{t,r) 



(81) 



= e-3(*'''\/,5 + L2£^^f^ + p2g2/(t,r) sinh^(i,r). 
V P^i,?-) 



(82) 



After introducing these expressions in (^6|) , (|7^) , a first order equation for ^ 
is easily obtained. 



^(t,r) = -e-2g(t,r) {Pj(i^j.)sinh^(t,r) +F^(i,r)cosh^(i,r)} 



(83) 



/ , / p-2f(t,r) 

F{t,r) = eS(*''^\/5 + L2^— + p2g2/(t,r)^ 

V P[t,r) 

which expands into the following expression for ^, 



(84) 



X {coshe [Sgr + L^^ (,g, - U - ^) + P^e^f (g, + /,) 

+ sinh e [<5gt + ^ (ft - /t - ^) + ^^'e'^ i9t + ft)]} , (85) 

where the dependence of the functions on t, r has been omitted in order to avoid 
a large equation. 

The equivalent equation for past-directed geodesies can be also obtained. 
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(86) 



Now we can check the equations to extract conditions for causal geodesies 
to be complete. No attention needs to be paid to the equations for z and </>, 
since they are just quadratures to be solved once t, r are known as functions of 
the afline parameter and cannot be singular unless the functions involved in the 
integration are singular and we have already imposed that the metric functions 
must be smooth. 

We have to prevent the coordinates t, r from tending to infinity at a finite 
value of the afline parameter. Also the coordinate r must not tend to zero at 
a finite value of the affine parameter when the constant of geodesic motion L 
is nonzero. When L is zero, the denominators that depend on p disappear and 
the expressions for t, f, ^ are not singular at r = 0. 

First of all we shall exclude the possibility of having arbitrarily large values 
of 1^1 . There are two cases to consider, 

• ^ > 0, ^ > : We shall impose a condition in order to prevent this 
situation from lasting too long by requiring that ^ changes sign for large 
values of the time coordinate t. Since we may have geodesies with any of 
the constants of motion equal to zero, the terms that multiply them in 
the numerator of ^ must be treated independently. For instance, for S, we 
need, 

< cosh^ + sinh^ ^f = sinh^ {gr + gt) + e^"' gr, (87) 

and similar conditions for L and P. These conditions can be fulfilled by 
requiring. 



gr+gt>0 

{g-f^\np)r + {g-f-\np)t>0 , (88) 
{g + f)r + {g + f)t > 

for large values of t and increasing r. And for the remaining terms that 
multiply e~^, we shall have further to impose that, 



gr 

{ {g-f-\np)r . (89) 

[ {g + I)r 

be either positive or at most of the same order as their respective terms 
in (|88| ) so that e^^ turns them negligible, also for large values of t and 
increasing r. 

With all these requirements ^ becomes negative before t may develop a 
singularity. 
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• ^<0, ^<0:Asin the previous case, we shall prevent ^ from growing 
unboundedly large. Since the unfriendly denominators depending on p, 
that are only dangerous in this case since r decreases, only appear when 
L 0, we shall just require for nonzero L, that ^ becomes positive for 
large values of t. Again, by imposing such restriction on the numerator of 
^ one finds, for large values of t and decreasing r, 



5{gt - 9r} + {{g-f- Inp)* -{g-f- Inp),.} + 

PV/{(5 + /),-(5 + /),,}>0 (90) 

by making use of the identity cosh^ = — sinh^ + e^. Since L 7^ 0, its term 
may help compensate for the others to achieve a positive numerator. The 
terms that multiply the exponential of ^, 



5gr (^gr - fr ^ j) + P'^'^ (dr + fr) , (91) 

must be either negative or at most of the same order as the term in ( pO[ ) 
for large values of t and decreasing r. 

Finally we shall prescribe a behaviour for t similar to the one used for lightlike 
radial geodesies to avoid faster than exponential growth for large values of t, 



git.r) 

g{t,r) + J{t,r) + \np ) > - \n\t + a\ + b, (92) 
g{t,r)- f{t,r) 

the first of which is obviously weaker than the condition already required for 
lightlike radial geodesies. The constants a, h need not be the same for the three 
inequalities. 

No similar restrictions need be imposed on r since this derivative cannot 
become singular before i does according to equation (|73|). 

Similar results can be derived for past-directed causal geodesies just replac- 



ing the sign of the derivatives with respect to t for the opposite in (88), (|90D, 
for small values of the time coordinate. The other conditions remain unchanged 
for small values of t since they where obtained for 

According to |^ , the spacetime is globally hyperbolic since from the deriva- 
tion of the results it follows that every lightlike geodesic intersects once and only 
once every hypersurface of constant time and therefore they are Cauchy hyper- 
surfaces. 

These results can be summarized in the following theorem. 

Theorem: A cylindrically symmetric diagonal metric in the form ( |64[ ) with 
metric functions /, g, p is future causally geodcsically complete if con- 
ditions (|79|), (H), (11), (|0l), (|l|), (H) are fulfilled for large values of 
t. 
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In addition, the spacetime is globally hyperbolic. 

A fact that is worthwhile mentioning is the relation for the orbit of a future- 
directed geodesic, 

dv 

— ==tanh^. (93) 

This means that, for events (to,ro), {t,r), tq < r, < on the geodesic we 
have, 

r<ro + it-to), (94) 

and if tq > r, to < t, then, 

r>ro-{t~to). (95) 

Furthermore, the orbit lies above its tangent lines while ^, ^ > and there- 
fore, we get a lower bound. For tq < r, to < t, 

r >ro+tanhC(to,ro) • (t-to), (96) 
And similarly, for ^, ^ < 0, rg > r, to < t, 

r < ro +tanhe(to,ro) • (t-to), (97) 

6 Examples 

In this section we shall show explicitly how the theorem works with all known 
diagonal nonsingular models. 

• Senovilla: This is the first known nonsingular cosmological model in 
the literature. It describes a universe in a radiation dominated epoch. 
Its geodesic completeness is proven in ||2^. The metric functions for this 
model are, 



g{t, r) = 2 Incosh(at) + lncosh(3ar), 
f{t, r) = lncosh(a t) + ^ lncosh(3 ar), 

p(t,r) = -!^cosh(ai)sinh(3ar)cosh"^/^(3ar), (98) 

where a is a positive constant related to the maximum of the pressure. 

Since they are even functions of t we need not worry about past-directed 
geodesies. 

All three functions in (|79|),(p2[) are positive except for constants and the 
term ln(sinh(3ar)) in the third function. For increasing r, it docs not 
mean a problem. For decreasing r, it is bounded by (Bq), 
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ln(sinh(3ar)) > ln(sinh(3a|i - ro + to\)), (99) 



and therefore they fulfil this condition. 



The first and third conditions in (88) are always satisfied for positive t 



and the same happens with the second for large values of r. 

The expressions in (^9|) do not depend on t and are always positive, except 
the second one, which requires large values of r. 

All three terms in (|9^) independently are positive for small values of r and 
positive t. 

Concerning (|9l]), the term as a whole is negative for small r and positive 
t. 

Hence this spacetime is geodesically complete. 

• Ruiz-Senovilla: [p3| This family includes the previous one as a subcase 
for K — I and n = 3. The metric is not written in isotropic form in the 
original paper, but it can be cast in that form by the change of variable, 



n a 



cosh(n a Ft) 



1 



P{x) = y/x^ + {K - l)a;(2"-i)/« - K, (100) 

where P{x) is a growing function for x > \. We shall need for our reasoning 
that, 

dR{r) ^ Pjcoshin a Rjr))) 
dr sinh(n a i?(r)) ' 

which tends to a positive constant when r tends to zero, to one when r 
tends to infinity and is always positive for positive r. 



1 Tl 7i — 1 

g{t,r) — — - — Incosh(ai) H — In cosh(n a i?(r)), 

TL — 1 Tl — 1 

f(t, r) = lncosh(a<) H In cosh(n a i?(r)), 

2 2 ri 

p(t,r) = cosh(ai) cosh(^""^/"(nai?(r))P(cosh^(nai?(r))), 

na L 

L = K-^. (102) 

The ranges of the parameters are n > 3, K > and a is again an arbitrary 
positive constant and every metric function is even in the coordinate time 
t. 
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Every function in in (|79|),(|92[) is either positive or involves a term that 
behaves as In r for small r and therefore can be bounded like in the previous 
example. 

The expressions in ( |8^ ) are all positive for large r and positive t as it 
happened in the Senovilla subcase. The first and third one are positive 
for all r. 

The radial derivatives in ( ^9|) are again positive for large values of r. The 
first and third one are positive for all r. 

Again the three terms in ( po| ) are positive for large t and small values of 



Finally, the angular momentum term in (91) is unboundedly negative for 
small r while the others are bounded and therefore the whole term is 
negative for small values of the radial coordinate. 

Hence the whole set of conditions is fulfilled and the family is geodesically 
complete. 

Mars: This is a family of nonsingular models with no equation of state 
and depending on two parameters a, related to the maximum density of 
perfect fluid, and c > 3/4,^ 



g{t,r) = -i lncosh(2at) + f cosh^(2 a i) sinh^(2 a r), 

f{t,r) — i lncosh(2ai) — Incosh(ar), 

1 sinh(2 ar) , , 

2 a cosh(2 a t) 

Again the metric functions are even in the time coordinate t. 

The negative terms in (^, ( p2| ) are counteracted by the exponential 
growth of the c terms for large r and t. 

Concerning (p8|), the three expressions are positive for positive t except 
for bounded terms in t and r. 

The same sort of behaviour is to be found in ( p9| ) but for all values of t. 

The terms in the expression of (|9^) are all bounded or tend to zero except 
for 2a coth(2ra). Therefore, this term is positive for small r. The reasoning 
why the terms on sinh(ar) tend to zero faster than the ones with cosh(at) 
is explained in the next example. 

Finally, the expression in (^TJ) is negative for small values of r since every 
term tends to zero or is boimded except for the — 2acoth(2ra) term. 

Therefore, this family is also geodesically complete. 



^Note that there is a missprint in the original paper. 
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• Fernandez- Jambrina: This non-separable metric, that has been obtained 
in section 4, corresponds to a spacetime filled with a stiff perfect fluid. 
Both parameters, a and /3, are positive, a is related to the maximum of 
the density of the fluid. The metric functions are, 



g{t,r) = r^ + 2Mr^ +pt^ + 2P^r^ e, 

/(t,r) = f(r2+2t2), p{t,r)^r, (104) 

and again they are even functions of time. 

The three functions in (|79|),(|9^) are positive, except for a constant and a 



term Inr to which applies the same reasoning used for Senovilla's metric 
and so this condition is trivially satisfied. 



The first and third expressions in (88) are positive for positive values of t. 
The second expression requires in addition large values of r. 

The same happens for (|8^), regardless this time of the sign of t. 

The terms in ( |90| ) are all either positive and growing linearly with t or as 
1/r or tending to zero. Note that the terms polynomic in r tend to zero 
although they may be multiplied by powers of t since according to (^7|) r 
decreases to zero faster than t grows to infinity. 

Finally, all of the terms in (|9^) tend to zero except which makes 

the expression negative for small values of r. 

Therefore these spacetimes are geodesically complete. 



7 Summary 

In this paper a new approach for G2 inhomogeneous cosmological models 
has been derived. The description is grounded on kinematical properties 
of spacelike congruences that have been encoded in several 1-forms. The 
exterior differential system, equivalent to Einstein equations, has been 
submitted to different kinds of simplifications. Based in this formalism 
several consequences are drawn. It is shown that when the transitivity 
surfaces are null a perfect fluid is not admissible. Also, when these surfaces 
are spacelike they are not trapped sets. 

As an application to non-singular cosmological models a solution is re- 
obtained within the formalism. A theorem is written stating a sufficient 
condition for diagonal models to be singularity-free and globally hyper- 
bolic. All known non-singular diagonal perfect fluid spacetimes fit in this 
framework. 

The present work has been supported by Direccion General de Ensenanza 
Superior Project PB95-0371. The authors wish to thank F. J. Chinea and 
M. J. Pareja for valuable discussions. 
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